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要 :【 目 的 】 本 研究 旨 在 通过 差异 表达 基因 (differentially expressed gene, DEG) 分 析 以 及 毒 力 因子 
和 其 他 侵 染 相关 因子 分 析 , 在 转录 组 水 平 揭示 东方 蜜蜂 微 孢 子 虫 Nosema ceranae 侵 染 意大利 蜜蜂 
Apis mellifera ligustica 的 分 子 机 制 。 【方法 】 基 于 前 期 已 获得 高 质量 的 东方 蜜蜂 微 孢 子 虫 纯化 孢子 
(NeCK) 及 侵 染 意大利 蜜蜂 工蜂 7 和 10 d HARA E SE COT & (PSA NecT1 和 NcT2) 转录 组 数 
据 , 根 据 P<0.05 H.llog, (Fold change) | z1 的 标准 ,通过 比较 分 析 筛 选 出 NecCK vs NeT1, NeCK vs 
NcT2 和 NcTI vs NcT2 比较 组 的 DEC。 通 过 相关 生物 信息 学 软件 对 上 述 DEG 进行 Venn 分 析 、GO 
分 类 和 KEGG 代谢 通路 富 te Cen at 

Oe Sb gm eg 过 RT-qPCR 验证 转录 组 数据 及 DEG 表达 

。【 结 果 ] 从 NeCK vs NcTI , NeCK vs NcT2 fe NcTI vs NcT2 比较 组 分 别 鉴定 出 1397, 1 497 和 

个 DEG。Venn 分 析 结 果 显 示 各 比较 组 共 A HL Alte TRAE BI 分 别 为 10 Fe 1 A, GO 分 类 结果 
显示 ,NcCK vs NcTI 和 NcCK vs NcT2 中 DEG 富 集 数 最 多 的 功能 条 目 为 代谢 进程 .细胞 进程 单 组 
织 进程 .细胞 、 细 胞 组 件 、 细 胞 器 、 众 化 活性 和 结合 ,而 NcTl vs NcT2 中 DEC 富 集 数 最 多 的 是 代谢 进 
程 、 细 胞 进程 . 单 组 织 进 程 、 众 化 活性 和 结合 。KEGG 代谢 通路 富 集 分 析 结 果 显 示 , NcCK vs NecTl 和 
NeCK vs NcT2 中 DEG 分 别 富 集 到 80 和 79 条 通路 ; 富 集 在 糖 酵 解 / 糖 异 生 和 MAPK 信号 通路 的 上 
调 基因 数量 多 于 下 调 基因 。 毒 力 因子 分 析 结 果 显 示 , 孢 壁 蛋 白 9 基因 和 和 孢 壁 蛋白 12 基因 在 NcCK 
vs NcTI 和 NcCK vs NcT2 中 均 下 调 表 达 , 孢 壁 蛋白 8 基因 仅 在 NecCK vs NeTI FRE TFH; 96 9 FA 

壁 蛋白 前 体 基 因 、 孢 壁 和 锚 定 盘 复合 蛋白 基因 、 几 丁 质 合 酶 基因 、 极 管 蛋白 基因 、 葛 麻 毒素 B 凝集 
素 基因 的 表达 水 平 在 NeCK vs NeT1 和 NeCK vs NeT2 中 表现 为 上 调 。 侵 染 相关 因子 分 析 结 果 表 明 ， 
糖 酵 解 途径 的 3 个 关键 酶 基因 在 NcCK vs NOTI 和 NcCK vs NeT2 中 上 调 表 达 ;3 个 涉及 ATP/ADP 
移 位 酶 的 基因 在 NcCK vs NcTl fe NcCK vs NcT2 中 上 调 表 达 , 但 有 1 个 表达 量 下 调 ;2 个 涉及 ABC 
转运 蛋白 的 基因 在 NcCK vs NcTl 和 NcCKw NcT2 中 上 调 表 达 ue TIT RT-qPCR 结果 
证 实 了 本 研究 中 转录 组 数据 及 DEG 表达 趋势 的 真实 可 靠 性 。 【结论 】 本 研究 通过 比较 分 析 解 析 东 
方 宫 蜂 微 孢 子 虫 侵 染 意大利 蜜蜂 工蜂 过 程 的 转录 组 动态 ,揭示 ieee. FE do LAZO 
白 、 几 丁 质 酶 、 极 管 蛋白 和 葛 麻 毒素 B 凝集 素 等 毒 力 因子 编码 基因 ,以 及 己 糖 激 酶 .丙酮 酸 激酶 6- 
磷酸 果糖 激酶 .ATPAADP 移 位 酶 和 ABC 转运 蛋白 等 侵 染 相关 因子 编码 基因 在 病原 增殖 中 扮演 重 
要 角色 ,为 阅 明 东方 蜜蜂 微 孢 子叶 的 侵 染 机 制 提供 了 基础 。 
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Transcriptomic analysis reveals the molecular mechanism underlying 


Nosema ceranae infection of Apis mellifera ligustica 

GENG Si-Hai*, ZHOU Ding-Ding*, FAN Xiao-Xue, JIANG Hai-Bin, ZHU Zhi-Wei, WANG Jie, FAN 
Yuan-Chan, WANG Xin-Rui, XIONG Cui-Ling, ZHENG Yan-Zhen, FU Zhong-Min, CHEN Da-Fu, 
GUO Rui* (College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry 
University, Fuzhou 350002, China) 

Abstract: [ Aim] This study aims to reveal the molecular mechanism underlying Nosema ceranae 
infection of Apis mellifera ligustica at the transcriptome level by investigation of differentially expressed 
genes ( DEGs) , virulence factors and infection-associated factors. [Methods] Following the criteria of 
P x:0.05 and llog, (Fold change) | =1, DEGs within NcCK vs NecT1 , NcCK vs NeT2 and NcTI vs NcT2 
comparison groups were screened out via the comparative analysis based on our previously obtained high- 
quality transcriptome datasets from purified spores of N. ceranae ( NCCK) and N. ceranae in the midgut 
of A. m. ligustica workers at 7 d and 10 d post infection ( NcT1 and NcT2, respectively). Venn 
analysis, GO classification and KEGG pathway enrichment analysis of DEGs were conducted using the 
related bioinformatic software, virulence factors and infection-associated factors of N. ceranae were 
summarized and analyzed based on Nr database annotations, KEGG database annotations and related 
literature documentations. The transcriptomic data and expression trends of DEGs were verified by RT- 
qPCR. [Results] In total, 1 397, 1 497 and 52 DEGs were identified in NcCK vs NeT1 , NcCK vs NcT2 
and NcT1 vs NcT2, respectively. Venn analysis showed that 10 up-regulated genes and one down- 
regulated gene were shared by various comparison groups. GO functional classification results showed the 
largest functional terms of DEGs in NcCK vs NcTI and NcCK vs NcT2 were metabolic process, cellular 
process, single-organism process, cell, cell part, organelle, catalytic activity and binding, while DEGs 
in NeT1 vs NcT2 were mostly enriched in metabolic process, cellular process, single-organism process, 
catalytic activity and binding. In addition, KEGG pathway enrichment analysis indicated that DEGs 
within NcCK vs NeT1 and NeCK vs NcT2 were enriched in 80 and 79 pathways, respectively. The 
number of up-regulated genes enriched in glycolysis/gluconeogenesis and MAPK signaling pathway was 
more than that of down-regulated genes. Investigation of virulence factors displayed that spore wall protein 
9 encoding gene and spore wall protein 12 encoding gene were down-regulated in NcCK vs NeT1 and 
NcCK vs NcT2, while spore wall protein 8 encoding gene was down-regulated only in NcCK vs NcT1. In 
additional, the expression levels of spore wall protein precursor encoding gene, spore wall and anchoring 
disk complex protein encoding gene, chitinase encoding gene, polar tube protein encoding gene, and 
ricin B lectin encoding gene were all up-regulated in NcCK vs NcTI and NcCK vs NcT2. Moreover, the 
infection-associated factor analysis demonstrated that in NcCK vs NcT1 and NcCK vs NeT2 three key 
enzyme genes engaged in glycolytic pathway were up-regulated, while in NcCK vs NcT1 and NeCK vs 
NcT2 three ATP/ADP transferase-associated genes were up-regulated but one down-regulated. Two ABC 
transporter-associated genes were up-regulated in NcCK vs NcT1 and NcCK vs NcT2, while four down- 
regulated. Finally, the RT-qPCR results verified the authenticity of the transcriptomic data and 
expression trends of DEGs. [Conclusion] In this study the transcriptomic dynamics of N. ceranae 
infecting the A. m.  ligustica worker was deciphered through the comparative analyses. Our findings 


revealed that genes encoding virulence factors including spore wall protein, spore wall and anchoring disk 
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complex protein, chitinase, polar tube protein and ricin B lectin, and genes encoding infection-associated 


factors such as hexokinase, pyruvate kinase, 6-phosphofructokinase, ATP/ADP transferase and ABC 


transporters, are likely to play key roles in pathogen proliferation, providing a basis for clarifying the 


infection mechanism of N. ceranae. 


Key words: Apis mellifera ligustica; Nosema ceranae; transcriptome; differentially expressed gene; 


infection mechanism; virulence factor; infection-associated factor 





i Hl FH ee Ph Sat E B AP JEL ad J 
iz ?$ ^E TEE oh 9] A Z6 T$ E 3 7] ( Cuomo et. al., 
2012), Á 1857 ERIE BRK AE Bombyx mori 中 
RIMET E Nosema bombycis ( Nageli, 1857) 
以 来 ,至 今 已 经 有 1 400 Z RE GUGCT EEE hK 
(吴杰 等 , 2017) 。 微 孢子 虫 的 基因 组 紧凑 ,缺少 典 
型 的 线粒体 ,取而代之 是 一 种 微粒 体 ,因而 高 度 依赖 
和 宿主 细胞 提供 能 量 ( Burri et al., 2006) 。 

ATT FE ME ELF E Nosema ceranae 是 专 性 寄生 
蜜蜂 中 肠 上 皮 细 胞 的 单 细 胞 真菌 病原 ,可 引起 蜜蜂 
微 孢 子 虫 病 ( Forsgren and Fries，2010 ) 。 Freis 等 
(1996) 首 次 在 北京 地 区 的 东方 蜜蜂 Apis cerana 蜂 
群 中 分 离 和 鉴定 出 东方 蜜蜂 微 移 子 虫 ,之 后 该 抱 子 
虫 逐渐 扩散 到 欧洲 ( Higes et al., 2006) 和 中 国 台 湾 
地 区 (Huang et al., 2007) 的 西方 蜜蜂 Apis mellifera 
峰 群 , 现 已 蔓延 至 世界 各 养 峰 国家 和 地 区 (Shumkova 
et al., 2018)。 前 人 的 研究 结果 表明 东方 蜜蜂 微 孢 
子 虫 感染 可 导致 成 年 蜜蜂 肠 道破 坏 TA A SCRA AL 
脂 质 代谢 加 快 和 脂 质 储存 减少 ,免疫 基因 和 细胞 调 
亡 被 抑制 ,能 量 胁迫 ,采集 日 龄 提前 及 采集 能 力 下 降 
等 ;东方 蜜蜂 微 孢 子 虫 感染 严重 可 造成 蜂 群 的 群 势 
降低 、 采 蜜 量 减少 以 及 越冬 失败 ( Higes et al., 2007; 
Antünez et al., 2009; Mayack and Naug, 2009; Botías 
et al., 2013; Chen et al., 2013; Mike et al., 2013; 
Kurze et al., 2015) 。 此 外 ,有 研究 表明 东方 密 蜂 微 孢 
子 虫 与 蜂 群 衣 演 失调 症 之 间 存 在 密切 的 关联 ( Higes 
et al., 2009; Simeunovic et al., 2014)。 东 方 蜜蜂 微 
孢子 虫 在 体外 以 休眠 态 的 孢子 存在 ,孢子 壁 由 几 丁 
质 和 孢 壁 蛋白 组 成 ,能 够 抵御 外 界 的 不 良 环境 ,孢子 
被 蜜蜂 摄取 之 后 ,在 蜜蜂 中 肠 的 特殊 理化 环境 的 刺 
激 下 发 芽 , 内 部 高 度 压缩 盘旋 的 极 丝 弹射 而 出 , 刺 穿 
中 肠 上 皮 细 胞 并 将 有 侵 染 性 的 胞 原 质 注入 ,随后 进 
入 繁殖 期 ,利用 宿主 的 物质 与 能 量 合成 系统 进行 增 
殖 , 孢 子 数量 逐渐 增多 ,最终 导 致 宿主 细胞 裂解 , 释 
放出 的 成 熟 孢 子 继续 感染 临近 的 上 皮 细 胞 ,或 者 随 凑 
便 排出 体外 ,感染 新 的 宿主 (Simeunovic et al., 2014)。 

随 着 高 通 量 测序 技术 和 生物 信息 学 算法 的 发 


























































































































展 ,RNA-Seq 技术 不 仪 广泛 应 用 于 人 类 ( Viollet et 
al., 2017) .动物 (李鹏 飞 等 ,2018 ) 和 植物 (Mun et 
al., 2017) 的 研究 ,在 微生物 的 相关 研究 中 也 得 到 了 
较 好 的 应 用 ,包括 绿 僵 菌 Metarhizium anisopliae 
(Wang et al., 2014) , Œ [C AX Trichoderma reesei 
(Ries et al., 2013) , Z E Git Hd F E ( Liu et al., 
2016a) All 2 EE ER SE Fe] Ascosphaera apis ( Sf R F, 
2017) 等。Huang (2016 ) XE Ie As 77 EE oh fd T 
虫 后 1 -6 d 的 西方 塞 蜂 中 肠 进 行 测序 ,从 混合 数据 
中 筛选 出 东方 蜜蜂 微 移 子 忠 的 数据 ,进而 通过 统计 
东方 蜜蜂 微 孢 子 虫 的 基因 表达 量 , 发 现 1 122 个 基 
因 可 分 为 4 种 表达 模式 ;随后 ,该 团队 对 东方 蜜蜂 微 
孢子 虫 的 Dicer 进行 RNAi, 然 后 通过 转录 组 测序 和 
分 析 发 现 Dicer 调控 东方 蜜蜂 微 孢 子 虫 的 基因 表达 
( Huang et al., 2018) 。 笔 者 团队 前 期 结合 链 特异 性 
建 库 方法 和 二 代 测 序 技术 对 东方 蜜蜂 微 抱 子 虫 孢 子 
进行 了 深度 测序 ,在 全 基因 组 水 平 预 测 、 鉴 定 和 分 析 
了 东方 蜜蜂 微 移 子 虫 的 长 链 非 编码 RNA (long non- 
coding RNA, IncRNA) 和 环 状 RNA (circular RNA, 
circRNA ) ( Guo et al., 2018a, 2018b) 。 但 对 于 东方 
AE EE ok PEL HY ZH OT FE BC BE NM US PR BON i I, 
DR Ta SENE NL EL T R AZAA GL tl BAS BR] A o 

FE A ARTH Cou] AR 77 S8 HRY 7 和 
10 d 的 意大利 蜜蜂 Apis mellifera ligustica 工蜂 中 肠 
进行 转录 组 深度 测序 ,通过 连续 比 对 西方 蜜蜂 和 东 
方 蜜蜂 微 孢 子 虫 的 基因 组 过 滤 得 到 东方 蜜蜂 微 孢 子 
虫 的 纯净 数据 ,结合 东方 蜜蜂 微 孢 子 虫 了 钨 子 的 转录 
组 数据 (Guo et al., 2018a) 开展 了 东方 蜜蜂 微 孢 子 
虫 的 高 表达 基因 (highly expressed gene, HEG ) 分 
析 , 解 析 了 东方 蜜蜂 微 孢 子 虫 侵 染 意大利 蜜蜂 过 程 
的 HEG 表达 谱 及 潜在 作用 ( 熊 染 玲 等 , 2019)。 本 
研究 在 前 期 研究 基础 上 对 东方 蜜蜂 微 人 雹 子 虫 进行 深 
入 的 转录 组 分 析 , 进 一 步 通过 分 析 病 原 的 毒 力 因子 
和 侵 染 相关 因子 ,揭示 东方 蜜蜂 微 孢 子 虫 侵 染 意 
利 蜜蜂 的 分 子 机 制 。 研 究 结果 可 为 探 明 东方 蜜蜂 微 
孢子 虫 的 侵 染 机 制 提供 基础 ,为 蜜蜂 微 雹 子 虫 病 的 
治疗 提供 潜在 靶 点 。 
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1 材料 与 方法 


1.1 转录 组 数据 来 源 

前 期 研究 中 ,笔者 团队 利用 基于 链 特 异性 建 库 
的 IncRNA-Seq 4 /ROVE AR 27 SE E idc HEY 2 (58 
子 进 行 了 深度 测序 (NeCK-1，NcCK-2 和 NeCK-3 为 
3 个 生物 学 重复 ), 获得 了 高 质量 的 mRNA 和 
IncRNA 组 学 数据 (Guo et al., 2018a) , 其 中 的 
mRNA 组 学 数据 可 作为 本 研究 中 的 对 照 组 数据 ;此 
外 ,还 利用 二 代 测 序 技术 对 东方 蜜蜂 微 伯 子 虫 感染 
7 d(NcTI-1, NeT1-2 和 NceT1-3 为 3 个 生物 学 重复 ) 
和 10 d(NeT2-1, NeT2-2 和 NeT2-3 293 个 生物 学 
复 ) 的 意大利 蜜蜂 工蜂 中 肠 进行 了 全 转录 组 测序 ， 
获得 了 高 质量 的 测序 数据 ( 熊 染 玲 等 , 2019) ,其 中 
的 mRNA 组 学 数据 包含 宿主 来 源 的 数据 和 病原 来 
源 的 数据 ,因此 通过 连续 比 对 核糖 体 数 据 库 .西方 蜜 
蜂 基因 组 和 东方 蜜蜂 微 移 子 忠 基因 组 对 上 述 混合 类 
据 进 行 过 滤 , 筛 滤 得 到 的 病原 的 mRNA 组 学 数据 
( 即 侵 染 意大利 蜜蜂 工蜂 的 东方 蜜蜂 微 移 子 虫 的 
mRNA 组 学 数据 ) 作为 本 研究 中 的 处 理 组 数据 。 

上 述 转录 组 测序 的 原始 数据 已 上 传 美国 国家 生 
物 技术 信息 中 心 (NCBI) SRA 数据 库 , Bioproject 号 
分 别 为 PRJNA395264 ( NcCK ) 和 PRJNA406998 
(NcTI 和 NeT2 ) 。 
1.2 差异 表达 基因 (DEG ) 分 析 

采用 FPKM 算法 计算 和 归 一 化 基因 表达 量 。 利 
用 edgeR 软件 筛选 NeCK vs NeT1, NeCK vs NeT2 和 
NcTl vs NcT2 比较 组 的 DEG, RME XW P <0.05 H 
llog, ( Fold change) | =1, Al) FA OmicShare 在 线 平 台 
(www. omicshare. com) 对 各 比较 组 中 上 调 基 因 和 下 
调 基因 进行 Venn 分 析 , 采 用 默认 参数 。 参 照 郭 窒 
等 (2017) 的 方法 ,利用 BLASTall 软件 将 DEG. 比 对 
到 GO 和 KEGG 数据 库 , 进 行 功 能 分 类 和 代谢 通路 
注释 ,采用 默认 参数 。 
13 ” 毒 力 因 子 和 侵 染 相关 因子 分 析 

AR Ty SEE LRL T B ee PE SEHE UAE SEE PES E 
皮 细 胞 的 真菌 病原 ,其 增殖 所 需 的 物质 和 能 量 完全 
依赖 宿主 细胞 提供 ( Burri et al., 2006)。 根 据 DEG 
在 Nr 数据库 的 注释 信息 ,同时 参考 前 人 关于 东方 蜜 
蜂 微 孢子 虫 的 基因 组 学 、 转 录 组 学 和 分 子 生 物 学 的 
人 研究 报道 ( Cornman et al., 2009; Paldi et al., 2010; 
Huang et al., 2016, 2018; Pelin et al., 2016; 
Rodríguez-García et al., 2018) WR Fh ih T HEY 

















































































































相关 研究 报道 ( Corradi and Slamovits, 2011; Katinka 
et al., 2011; Nakjang et al., 2013; Simeunovic et al., 
2014) fith JL T EAE He Re SR A, BE REE BORE 
集 素 、 极 管 蛋 白 等 毒 力 因子 相关 DEG, 以 及 己 糖 激 
酶 .丙酮 酸 激 酶 .磷酸 果糖 激酶 .ATP/ADP 移 位 酶 、 
ABC 转运 蛋白 等 侵 染 相关 因子 相关 DEG ,并 统计 上 
述 相 关 DEG 在 各 组 样品 中 的 表达 量 (FPKM ) 及 在 
各 比较 组 中 的 log, (Fold change) 。 
1.4 DEG 的 RT-qPCR 验证 

从 NeCK vs NcT1, NeCK vs NcT2 和 NcT1 os 
NeT2 比较 组 中 分 别 随机 选取 7 个 DEC ,其 中 NeCK 
vs NeTI. 中 包含 6 个 上 调 基因 和 1 个 下 调 基 因 ， 
NcCK vs NcT2 中 包含 6 个 上 调 基 因 和 1 个 下 调 基 
[A ,NeTI vs NeT2 中 包含 5 个 上 调 基因 和 2 个 下 调 
基因 ,以 actin (GenBank 登录 号 : 36320842 ) 作为 内 
参 基因 ,利用 Primer Premier 5 设计 引物 ( 表 1) , 委 
托福 州 擎 科 生 物 有 限 公 司 合成 。 利 用 RNA 提取 试 
WA (TaKaRa 公司 , 大连) 分 别提 取 NeCK 以 及 
NcTI 和 NeT2 中 肠 总 RNA ( RERI 4$, 2019) ,然后 
利用 cDNA 第 1 链 合成 试剂 盒 ( Vazyme 公司 ,南京 ) 
反 转 录 成 cDNA, 作为 模板 进行 qPCR, 反 应 在 
QuantStudio 3 ( ThermoFisher 公司 ,美国 ) 上 进行 。 
qPCR 反应 按照 SYBR Green Dye 试剂 盒 ( Vazyme 公 
司 ,南京 ) 操 作 说 明 书 进行 ,本 实验 进行 3 次 生物 学 
重复 ,每 个 生物 学 重复 包含 6 头 意大利 蜜蜂 工蜂 中 
肠 。 反 应 体系 (20 wh) : 正 反 向 引物 (2.5 pmol/L) 
各 1 uL, cDNA 模板 1 uL, SYBR Green Dye 10 uL, 
DEPC 水 7 kwL。 反 应 程序 : 95°C 预 变性 min; 95°C 
变性 15 s，60Y 延伸 30 s, 共 40 个 循环 。 通 过 
2 “全 法 对 上 述 基 因 的 相对 表达 量 进 行 计 算 。 利 用 
Graph Prism 7 软件 进行 数据 分 析 和 绘图 。 
































2 结 


2.1 东方 蜜蜂 微 孢 子 虫 转录 组 的 DEG 分 析 

差异 分 析 结 果 显 示 , NeCK vs NeT1, NeCK vs 
NcT2 和 NeT1 vs NcT2 比较 组 分 别 得 到 1 397, 1 497 
和 52 个 DEG ,其 中 上 调和 下 调 基因 分 别 为 497 和 
900 个 , 517 和 977 个 ,38 和 13 M, Venn 分 析 结 果 
显示 ,有 10 个 基因 在 上 述 3 个 比较 组 中 共同 上 调 表 
达 ( 图 1: A), 它 们 涉及 碳水 化 合 物 代谢 、 节 麻 毒素 
B 凝集 素 ,细胞 表面 糖 基 化 蛋白 \ 信 号 肽 酶 和 假定 蛋 
白 ( 表 2); 仅 有 1 个 下 调 表 达 的 假定 蛋白 编码 基因 
为 各 比较 组 所 共有 (图 1: B; 表 2)。 
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表 1 本 研究 使 用 的 引物 
Table 1 Primers used in this study 


GenBank 登录 号 


GenBank accession no. 


引物 序列 (5' -3") 产物 大 小 (bp) 


Primer sequence Product size 


Nr 注释 


Nr description 






































big F: ATGAAGACATTGGTGGGTTA m 26S 和 蛋白酶 调节 亚 基 6a 
R; CTGTTCCAGGAGGACCATAC 26S protease regulatory subunit 6a 
F: AAACACCAAAGAAACGAGTA 假定 蛋白 
9422527 101 
R; ATCTGATTCATCTTTTACCC Hypothetical protein 
F: CACCTGGAAGCAAACACTAA 交 运 动 蛋 F 
9424439 145 moo 
R: AAGAAGCACGCCACTCAAAA Nuclear movement protein 
F: TAACACGTAAAGCCAGCAAA : 亚 基 -7 
9422132 157 EAR p Ta 
R: CCAGTATCGTCAACACCTCC Proteasome subunit beta type-7 
F: CAGACTGTTGGTGCCGATGA 组 蛋白 h2b 
9424049 130 
R: ATITAGAAGGCTTAGATGGT Histone h2b 
9423882. F: TATGAAGCGTCCAATACCAA 3o CCAAT 结合 转录 因子 o. 亚 基 
R; CAAGTTCAGCATAAACACCC CCAAT binding transcription factor subunit a 
9422528 F: CTGCAATTATGTTGGGTAAA (55 Jte HSE A 
R; TATCTTGGTCTAAAGCGTGT Adenylate kinase related protein 
F: TAGTTAAGGGATTTGGTATG EMBER B 凝集 素 
9422645 176 人 
R: TATTAGTGATTCGTTATTCG Ricin B lectin 
F: TGCATAAGGTCTCGGGAAAT 假定 蛋白 
9423520 121 
R; TTACTCCATTGTGGGCGAAG Hypothetical protein 
F: TCCACCACTGCTAAACCATC 4 WEE 
9423191 : 135 TERIERA 
R: AGATTTCTCGGTAAGCTCCT Cell surface glycoprotein 
9424131 F: GCCTTGTGCCAACGGTAAAT iu 谷 氨 酸 - 半 胱 氨 酸 连接 酶 催化 亚 基 亚 型 1 
R; TCAAATGGAGAACACCTGCT Glutamate-cysteine ligase catalytic subunit isoform 1 
F: ACGATGTGGAGGTAGGAGAA XE E E 
9424136 100 Fon 
R; AGGCTGGCAGAATGTAATGT Hypothetical protein 
F: AAGTCAAGTCGGTCAAAGAG XE AF 
9422787 : 195 假定 蛋白 
R: CATGGATCGGATTCATAAGC Hypothetical protein 
0 F: AACTAGGTGCTGTTGGTAAA ibs 分 裂 和 聚 腺 彰化 因子 (CF-IA) 组 分 Rnal5 
R: TTCTGCTATGAGGGGTGATT Cleavage and polyadenylation factor ( CF-IA) component Rnal5 
F: CAAAGCCTGTGATGTGGA ON aE E 
36320842 A 106 MIRA 
R: GCCAAGTTTCTGATTTCCC Actin 
A B 
NcCK vs NcT2 
NeCK vs NcTI NcCK vs NcT2 


Fig. 1 





NeT1 vs NcT2 






NeCK vs NcTI 


NeT1 vs NeT2 








图 1 








各 比较 组 中 东方 密 蜂 微 爷 子 虫 转录 组 差异 表达 基因 (DECs) Venn 分 析 


Venn analysis of differentially expressed genes ( DEGs) in Nosema ceranae transcriptome between various comparison groups 


A: 上 调 基 因 Up-regulated genes; B: 下 调 基因 Down-regulated genes. NcCK ; 微 孢 子 虫 纯 化 的 孢子 Purified spores of N. ceranae; NcT1 , NcT2: 分 别 


为 感染 意大利 蜜蜂 工蜂 7 和 10 d 中 肠 





infection, respectively. 下 同 The same below. 





PY AY ZR Zr 88 E i fL T H 


E N. ceranae in the midgut of Apis mellifera ligustica workers at 7 and 10 d post 
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表 2 各 比较 组 中 东方 蜜蜂 微 孢 子 虫 转录 组 共同 上 调 与 下 调 基 因 信 息 统计 
Table2 Summary of commonly up- and down-regulated genes in Nosema ceranae transcriptome 
between various comparison groups 
GenBank NcCK 中 NeT1 中 NeT2 中 NcCK vs NcTI NcCK vs NcT2 NeT1 vs NcT2 
登录 号 FPKM 值 ”FPKM 值 FPKM fü log (差异 倍数 ) log, (差异 倍数 ) ”logz (差异 倍数 ) Nr 描述 
CenBank FPKM value FPKM value FPKM in NcCK vs NcTI NcCK vs NcT2 NcT1 vs NcT2 Nr description 
accession no. in NcCK in NcTI NcT2 log; (Fold change) log, (Fold change) log, ( Fold change) 
, 糖 基 化 蛋 F 
9423191 64 213 532 1.74 3.06 1.32 由 胞 表面 糖 基 化 重 白 
Cell surface glycoprotein 
PE EL ge : 
9424632 18 174 355 3.29 4.33 1.03 KERERE , 
Signal peptidase-like protein 
丙酮 酸 脱 氧 酶 el 组 分 B 亚 基 
9423243 59 342 869 2.53 3.88 1.34 Pyruvate dehydrogenase el 
component beta subunit 
9423279 28 163 328 2.52 3.54 1.01 磷酸 甘油 酸 激酶 
Phosphoglycerate kinase 
9423592 13 184 397 3.7 4.90 1.11 ERREUR 
Sugar kinase 
甘油 3- 磷酸 脱氧 酶 
9423035 72 431 973 2.58 3.76 1.17 Glycerol-3-phosphate 
dehydrogenase 
EREA B 凝集 素 
9422649 98 2 655 6 200 4.76 5.98 1.22 ERRER B MERR 
Ricin B lectin 
EREA B 凝集 素 
9422645 95 1 024 2 546 3.42 4.74 1.31 ERRER B MERR 
Ricin B lectin 
EREA B 凝集 素 
9422719 30 901 2 437 4.88 6.32 1.46 ERRER B MERR 
Ricin B lectin 
9422567 23 1 883 3 990 6.37 7.45 1.08 POESIE 
Hypothetical protein 
9422077 411 68 26 -2.59 -3.99 -1.39 POESIE 
Hypothetical protein 
2.2 东方 蜜蜂 微 孢子 虫 转录 组 中 DEG 的 GO 分 进一步 对 各 比较 组 的 DEG 进行 KEGG 通路 富 


类 及 KEGG 通路 富 集 分 析 

GO 分 类 结果 表明 ,NcCK vs NeTI 中 DEG 涉及 
生物 学 进程 .细胞 组 分 和 分 子 功能 相关 的 23 个 功能 
条 日 ,其 中 代谢 进程 (221) 催化 活性 (195) .结合 
(193 )、 细 胞 进程 (193)、 细 胞 (116)、 细 胞 组 件 
(116) . 单 组 织 进程 (97 ) 和 细胞 器 (63 ) 的 条 目 富 集 
基因 数 最 多 ; NcCK vs NeT2 中 DEG 也 涉及 23 个 功 
能 条 目 , 同样 有 较 多 的 DEG 富 集 在 代谢 进程 
(218 )、 催 化 活性 (196)、 细 胞 进程 (195 )、 结 合 
(193) 细胞 (119)、 细 胞 组 件 (119)、 单 组 织 进程 
(94) 和 细胞 器 (51) 等 ;NecT1 vs NcT2 中 DEG 涉及 
生物 学 进程 (4) 、 细 胞 组 分 (5) 和 分 子 功 能 (2) 相关 
的 11 个 功能 条 目 , 其 中 催化 活性 、 细 胞 进程 .代谢 进 
TE . 单 组 织 进 程 和 结合 的 条 目 富 集 基 因数 最 多 ,分别 
为 11, 10, 10,7 和 6 个 (图 2)。 






































集 分 析 ,结果 显示 NeCK vs NeT1 中 有 355 个 DEG 
富 集 在 80 条 通路 ,涉及 代谢 .遗传 信息 加 工 .环境 信 
息 加 工 、 细 胞 进程 和 有 机 体系 统 等 5 个 大 类 ,其 中 富 
集 基 因数 最 多 的 是 新 陈 代 谢 途 径 (93) .次 级 代谢 产 
物 合成 (37) .核糖 体 (33) 抗生素 生物 合成 (29 ) 和 
真 核 生 物 核糖 体 生物 合成 (27) (R 3); NeCK vs 
NcT2 中 有 356 个 DEG 富 集 在 79 条 通路 ,其 中 富 集 
基因 数 最 多 的 是 新 陈 代谢 途径 次 级 代谢 产物 合成 、 
核糖 体 .抗生素 合成 和 真 核 生 物 核糖 体 生 物 合成 ,分 
BAY 92, 39, 32, 31 和 26 个 ( 表 3)。 此 外 NcCK vs 
NcT1 中 有 7 个 DEG 富 集 在 MAPK 信号 通路 ,其 中 
有 5 个 基因 上 调 表达 , GenBank. 登录 号 分 别 为 
9423510, 9424462, 9422987 ，9423963 和 9422306 ; 
另 有 2 个 基因 下 调 表 达 , GenBank 登录 号 分 别 为 
9422598 和 9424007 , NcCK vs NcT2 中 同样 也 有 7 个 
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NcCK vs NcTI 

转运 活性 Transporter activity 
结合 Binding 

代谢 进程 Metabolic process 
催化 活性 Catalytic activity 


细胞 器 组 件 Organelle part 

细胞 膜 组 件 Membrane part 

细胞 膜 Membrane 

大 分 子 复 合 物 Macromolecular complex 
细胞 器 Organelle 

细胞 Cell 





外 "| 
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NcCK vs NcT2 
转运 活性 Transporter activity 
结合 Binding | 
催化 活性 Catalytic activity 
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细胞 器 组 件 Organelle part — 


细胞 膜 组 件 Memb: art SSF 
细胞 肛 Membrane n 7 
大 分 子 复合 物 Macromolecular complex 物 
细胞 器 Organelle N 
细胞 Cell NN 4 


NcT1 vs NcT2 
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结合 Binding 
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催化 活性 Catalytic activity 











上 调 基因 Up-regulated gene 

调 基 因 Down-regulated gene 
生物 学 进程 Biological process 
分 子 功能 Molecular function 
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生物 学 调控 Biological regulation 
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细胞 成 分 组 成 或 生物 合成 Cellular component organization or biogenesis 
WUZ Response to stimulus 

细胞 组 件 Cell part 





代谢 进程 Metabolic process 


v 


包 进 程 Cellular process 







组 织 进程 Single-organism process 
学 调控 Biological regulation 
位 Localization 
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细胞 组 件 Cell part 
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生物 学 调控 Biological regulation 


图 2 各 比较 组 中 东方 蜜蜂 微 雹 子 虫 转录 组 差异 表达 基因 (DECs) 的 GO 分 类 


Fig. 2 GO categories of differentially expressed genes ( DEGs) in Nosema ceranae transcriptome 


between various comparison groups 


DEG 富 集 在 MAPK 信和 号 通路 ,包括 5 个 上 调 基 因 ， 
GenBank 登录 号 分 别 为 9423510, 9424610, 
9424462 , 9422987 和 9422306; 以 及 2 个 下 调 基 因 ， 
GenBank 登录 号 分 别 为 9411598 和 9424007。 在 
NcCK, NcTI 和 NcT2 中 , GenBank 登录 号 为 
9424462 的 基因 的 FPKM 值 分 别 为 966, 829 8122, 
GenBank 登录 号 为 9422306 的 基因 的 FPKM 值 分 别 
为 209, 205 和 5。 
2.3 ”东方 蜜蜂 微 孢 子 虫 的 毒 力 因子 分 析 

根据 NeCK vs NeT1, NeCK vs NcT2 和 NeT1 vs 
NeT2 的 DEG 的 Nr 数据 库 注释 情况 对 东方 蜜蜂 微 
孢子 虫 的 毒 力 因子 进行 统计 ,分 别 发 现 5 AP LEER 
白 基因 1 个 孢 壁 蛋白 前 体 基 因 和 1 4S TELE RIVER XE 
AE AE AEA , HPA 9 基因 ( GenBank 登 
录 号 : 9422782 ) Æ NcCK vs NcTI 和 NcCK vs NeT2 





中 下 调 表达 [log: (Fold change) 值 分 别 为 -2. 10 和 
-2.06],48 EE Æ A 12 4 DX] ( GenBank 登录 号 : 
9422437 ) 也 下 调 表达 [log, ( Fold change ) 值 分 别 为 
-0.86 FI -0.67 | ; EEE H 8 基因 (GenBank 登录 
号 : 9422311) E NcCK vs NeT1 中 表达 量 下 调 [log, 
(Fold change) 值 为 2. 05 ] ,但 在 NeCK vs NcT2 中 
其 表达 量 不 变 ; 在 NcCK vs NcTI 和 NcCK vs NcT2 中 
均 上 调 表 达 的 基因 包括 2 A fü EE 4 F1 3E DS 
[ GenBank 登录 号 为 9422653 的 基因 的 log, (Fold 
change ) 值 分 别 为 7.11 和 7.12; GenBank 登录 号 为 
9423533 的 基因 的 log, ( Fold change) 值 分 别 为 6. 32 
All 6.93 | 1 个 抱 壁 蛋白 前 体 基因 [ GenBank 登录 号 : 
9422389 log, ( Fold change) 值 分 别 为 5.03 和 5.46] 
和 1 个 孢 壁 和 锚 定 盘 复 合 和 蛋白 基因 | GenBank 登录 
号 :9424403,, log, (Fold change) 值 分 别 为 5.54 和 
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表 3 NcCK vs NcTI 和 NeCK vs NcT2 中 差异 表达 基因 ( DEGs) 富 集 数 前 18 的 KEGG 通路 
Table3 Top 18 KEGG pathways enriched by differentially expressed genes (DEGs) in NcCK vs NcT1 and NcCK vs NcT2 









































通路 NeCK vs NcTI NeCK vs NcT2 

Pathway 上 调 Up-regulated ”下 调 Down-regulated EJ] Up-regulated FY] Down-regulated 
新 陈 代 谢 途 径 Metabolic pathway 46 47 45 47 
次 生 代 谢 产物 的 生物 合成 
Biosynthesis of secondary metabolites A w p T 
核糖 体 Ribosome 16 17 14 18 
抗生素 的 生物 合成 Biosynthesis of antibiotics 17 12 20 11 
真 核 生物 核糖 体 的 生物 发 生 
Ribosome biogenesis in eukaryotes $ 5 à 和 
细胞 周期 -酵母 Cell cycle-yeast 10 16 8 17 
IC Hf Purine metabolism 8 13 7 14 
RNA 转运 RNA transport 5 16 4 17 
喀 喧 代谢 Pyrimidine metabolism 10 10 8 11 
mRNA 监测 通路 mRNA surveillance pathway 6 14 6 13 
减 数 分 型 -酵母 Meiosis-yeast 8 11 7 14 
内 质 网 蛋白 加 工 
Protein processing in endoplasmic reticulum x g = 
泛 素 介 导 的 蛋白 水 解 Ubiquitin-mediated proteolysis 4 13 4 13 
碳 代谢 Carbon metabolism 10 6 13 5 
蛋白 酶 体 Proteasome 13 2 13 4 
PET WR VIVRE Nucleotide excision repair 5 9 3 iô 
剪接 体 Spliceosome 3 10 4 10 
RNA 降解 RNA degradation 5 8 6 8 


4. 61]。 此 外 , 男 有 1 个 内 切 儿 丁 质 酶 基因 
[ GenBank 登录 号 : 9422749 , log, ( Fold change ) 值 分 
别 为 2. 94 和 3. 60]、1 个 几 丁 质 合 酶 1 基因 
[ GenBank 登录 号 为 9422288 的 基因 的 log ( Fold 
change) 值 分 别 为 6.08 和 5.56] 3 个 极 管 蛋 白 基 因 
[ GenBank 登录 号 为 9423198 的 基因 的 log, ( Fold 
change ) 值 分 别 为 5.62 和 6.34; GenBank 登录 号 为 
9423201 的 基因 的 log, (Fold change) 值 分 别 为 5.74 
和 6.37; GenBank 登录 号 为 9422340 的 基因 的 log, 
(Fold change ) 值 分 别 为 3.94 和 4.39] 和 7 NER 
毒素 B 凝集 素 基 因 | GenBank 登录 号 为 9422649 的 
基因 的 log, (Fold change) 值 分 别 为 4.76 和 5. 98; 
GenBank 登录 号 为 9422719 的 基因 的 log, ( Fold 
change ) 值 分 别 为 4.89 和 6.32; GenBank 登录 号 为 
9422642 的 基因 的 log, (Fold change ) 值 分 别 为 4.77 
和 4.98; GenBank 登录 号 为 9423906 的 基因 的 log, 
(Fold change) 值 分 别 为 2.32 和 3. 23; GenBank 登 
录 号 为 9422645 的 基因 的 log, ( Fold change) 值 分 别 
为 3.43 和 4.74; GenBank 登录 号 为 9422644 的 基 
因 的 log, (Fold change ) 值 分 别 为 2.72 和 3.21; 








GenBank 登录 号 为 9422647 的 基因 的 log, ( Fold 
change ) 值 分 别 为 5.62 和 3.38] 在 NcCK vs NcTI 和 
NcCK vs NcT2 中 均 上 调 表 达 ( 表 4)。 
2.4 ”东方 蜜蜂 微 孢 子 虫 的 侵 染 相关 因子 分 析 
进一步 对 东方 蜜蜂 微 和 抱 子 虫 的 其 他 侵 染 相关 因 
子 进行 统计 ,发 现在 NcCK os NcT1 和 NcCK vs NcT2 
中 均 上 调 表 达 的 糖 醉 解 关键 酶 基因 包括 己 糖 激酶 基 
[ GenBank 登录 号 : 9423428 , log, (Fold change) fH 
分 别 为 3. 38 和 2. 79 ] .丙酮 酸 激酶 基因 [ GenBank 
登录 号 : 9424331, log, (Fold change) [H7 32g 1. 81 
fll 2.69 | fll 6- ERRARE AA SE A [ GenBank 登录 号 : 
9424162, log, (Fold change) 值 分 别 为 1.91 $12.79] 
等 ( 表 5)。 TE NcCK vs NcTI 和 NcCK vs NcT2 中 均 
上 调 表达 的 有 3 个 ATP/ADP 移 位 酶 的 编码 基 
[ GenBank 登录 号 为 9424586 的 基因 的 log, ( Fold 
change) 值 分 别 为 2.31 和 2.34; GenBank 登录 号 为 
9422880 的 基因 的 log, (Fold change) 值 分 别 为 3. 60 
和 3.44; GenBank 登录 号 为 9424363 的 基因 的 log, 
(Fold change) 值 分 别 为 0.72 和 0.45]。 此 外 ,有 
1 4SATP/ ADP 移 位 酶 的 编码 基因 [| GenBank 登录 号 : 
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表 4 东方 蜜蜂 微 孢 子 虫 的 毒 力 因子 统计 
Table 4 Summary of virulence factors of Nosema ceranae 

GenBank NcCK 中 NcT1 中 NcT2 中 NcCK vs NcTI NcCK vs NcT2 

登录 号 FPKM 值 FPKM 值 FPKM 值 logs( 差 异 倍数 ) logs (差异 倍数 ) Nr 描述 

GenBank FPKM of FPKMof FPKM of NeCK vs NeT1 NcCK vs NcT2 Nr description 

accession no. NeCK NeTI NeT2 log, (Fold change) log, (Fold change) 

9422749 21 160 253 2.94 3.60 内 切 几 丁 质 酶 Endochitinase 

9422288 10 709 493 6.08 5.56 几 丁 质 合 酶 1 Chitin synthase 1 

9422389 170 5 575 7 481 5.03 5.46 孢 壁 蛋白 前 体 Spore wall protein precursor 
9422782 12 059 2 807 2 885 -2.10 -2.06 JE ETE 9 Spore wall protein 9 

9422311 2 797 675 2 794 -2.05 0.00 fü EAE |] 8 Spore wall protein 8 

9422431 1 618 894 1 018 - 0. 86 -0.67 $f REAE 12 Spore wall protein 12 
9422653 139 19 211 19 453 7.11 7.12 $f EAE H Spore wall protein 

9423533 131 10 480 16 021 6.32 6.93 TUEZ H Spore wall protein 

9424403 5 241 127 5.54 4.61 有 

Spore wall and anchoring disk complex protein 

9422649 98 2 655 6 199 4.76 5.98 HIRE Z B 凝集 素 Ricin B lectin 
9422719 30 901 2 437 4. 89 6.32 HIREK B 凝集 素 Ricin B lectin 
9422642 25 682 790 4.77 4.98 HIREK B REGEX Ricin B lectin 
9423906 260 1 297 2 438 2.32 3.23 BURG SE B 凝集 素 Ricin B lectin 
9422645 95 1 024 2 546 3.43 4.74 EERIE ZS B 凝集 素 Ricin B lectin 
9422644 16 103 145 2.72 3.21 EERIE ZS B 凝集 素 Ricin B lectin 
9422647 1 37 8 5.62 3.38 HIRE Z B 凝集 素 Ricin B lectin 
9423201 335 17 919 27 591 5.74 6.37 极 管 蛋白 2 Polar tube protein 2 

9423198 171 8 407 13 854 5.62 6.34 极 管 蛋白 1 Polar tube protein 1 

9422340 190 2 923 3 994 3.94 4.39 极 管 蛋白 Polar tube protein 

表 5 东方 窄 蜂 微 孢 子 虫 的 侵 染 相关 因子 统计 
Table 5 Summary of infection-associated factors of Nosema ceranae 

GenBank NcCK 中 Ne TI 中 Ne T2 中 NcCK vs NcTI NcCK vs NcT2 

登录 号 FPKM 值 FPKM 值 FPKM 值 logs( 差 异 倍数 ) logy (差异 倍数 ) Nr 描述 

GenBank FPKM in FPKMin FPKM in NeCK vs NcTI NeCK vs NeT2 Nr description 

accession no. NeCK NcTI NeT2 log, (Fold change) log, ( Fold change) 

9423428 14 144 96 3.38 2.79 己 糖 激酶 Hexokinase 

9424331 119 417 768 1.81 2.69 丙酮 酸 激酶 Pyruvate kinase 

9424162 152 572 1 051 1.91 2.79 6- 磷 酸 果糖 激酶 6-Phosphofructokinase 
9424586 768 3 818 3 885 2.31 2.34 ATP/ADP 移 位 酶 ATP/ ADP translocase 
9422880 86 1 039 928 3.60 3.44 ATP/ADP 移 位 酶 ATP/ ADP translocase 
9422375 8 212 655 551 -3.65 -3.90 ATP/ADP 移 位 酶 ATP/ADP translocase 
9424363 82 134 112 0.72 0. 45 ATP/ADP 移 位 酶 ATP/ ADP translocase 
9423050 5 764 458 568 -3.65 -3.34 ABC 转运 蛋白 ABC transporter 
9423041 62 221 224 1.83 1.85 ABC 转运 蛋白 ABC transporter 
9422518 60 115 94 0.93 0. 64 ABC 转运 蛋白 ABC transporter 
9423588 390 118 357 -1.72 -0.13 ABC 转运 蛋白 ABC transporter 
9422298 281 75 95 -1.91 -1.56 ABC 转运 蛋白 ABC transporter 
9422297 329 297 258 -0.19 -0.35 ABC 转运 蛋白 ABCtransporter 

9423588 83 48 132 -0.79 0.67 ABC 转运 蛋白 ABC transporter 





























3 期 耿 四 海 等 : 转录 组 分 析 揭 示 东 方 蜜蜂 微 抱 子 虫 侵 染 意大利 蜜蜂 的 分 子 机 制 303 











9422375, log, (Fold change) 值 分 别 为 - 3. 65 和 
-3.90] 在 NcCK vs NcTI 和 NcCK vs NeT2 中 均 下 
调 表 达 。 男 发 现 7 个 DEG 涉及 ABC 转运 蛋白 ,其 
中 有 2 个 | GenBank 登录 号 为 9423041 的 基因 的 
log, (Fold change) (H7 3] 2j 1. 83 和 1.85; GenBank 
登录 号 为 9422518 的 基因 的 log, ( Fold change) 值 分 
别 为 0.93 Al 0. 64 ] Æ NcCK vs NcTI 和 NcCK vs 
NcT2 中 均 上 调 表达 ;还 有 1 个 DEG( GenBank 登录 
号 : 9423588) 在 NeCK vs NeT1 中 下 调 表达 [log， 
(Fold change) 值 为 -0.79] ,而 在 NcCK vs NcT2 中 
表达 量 上 调 [log, (Fold change) 值 为 0.67] ;其 余 的 
4 个 ABC 转运 蛋白 基因 |[ GenBank 登录 号 为 
9423050 的 基因 的 log, ( Fold change ) 值 分 别 为 
-3.65 和 -3.34; GenBank 登录 号 为 9423588 的 基 
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图 3 KERHA P ARZT BE ME il PH RAE EIA ETA ( DEGs) 的 RT-qPCR 验证 


因 的 log, (Fold change) 值 分 别 为 -1.72 $1 - 0.13; 
GenBank 登录 号 为 9422298 的 基因 的 log, ( Fold 
change ) 值 分 别 为 -1.91 和 -1.56; GenBank 登录 
号 为 9422297 的 基因 的 log (Fold change) 值 分 别 为 
-0.19 FI -0.35 | E NcCK vs NcTI 和 NcCK vs NcT2 
中 表现 为 表达 量 下 调 ( 表 5)。 
2.5 ”东方 蜜蜂 微 孢子 虫 转录 组 DEG 的 RT-qPCR 
为 验证 RNA-Seq 结果 ,在 NeCK vs NecT1, NeCK 
vs NcT2 和 NcT1 vs NeT2 比较 组 中 分 别 随机 挑 取 7 
个 DEG 进行 RT-qPCR 验证 ,实验 结果 与 转录 组 数 
据 的 差异 表达 趋势 一 致 (图 3) ,证 明了 测序 数据 和 
基因 差异 表达 情况 的 可 靠 性 。 
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Fig. 3 


RT-qPCR verification of differentially expressed genes ( DEGs) in Nosema ceranae transcriptome 


between various comparison groups 
A: NeCK vs NeT1; B; NcCK vs NcT2; C; NeTl vs NeT2. 误差 棒 表 示 各 DEG 的 RT-qPCR 结果 的 方差 。Error bars represent the variance of RT- 


qPCR results of each DEG. 


3 讨论 


前 人 研究 报道 东方 蜜蜂 微 抱 子 忠 在 感染 西方 蜜 
蜂 后 10 -20 d 内 都 处 于 一 个 持续 增殖 的 过 程 ,在 此 
过 程 病原 的 孢子 数 持 续 增 多 ( Huang and Solter, 














2013 ) 。 前 期 研究 中 ,我 们 发 现 东 方 蜜蜂 微 孢 子 虫 
感染 的 意大利 蜜蜂 工蜂 的 死亡 率 在 感染 后 1 - 10 d 
范围 内 不 断 上 升 , 且 与 未 感染 工蜂 的 死亡 率 差 异 极 
显著 (Chen et al., 2019) ,表明 随 着 宿主 中 肠 内 东方 
蜜蜂 微 孢 子 虫 的 不 断 增殖 ,病原 数量 持续 增多 ,对 和 宿 



































304 昆虫 学 报 Acta Entomologica Sinica 63 3$ 





主 的 胁迫 压力 持续 上 升 。 在 另 一 项 研究 中 ,我 们 发 
现 东方 蜜蜂 微 孢 子 虫 的 孢子 数 在 4 -13 d 的 范围 内 
持续 增加 (未 发 表 数据 ) 。 本 研究 是 属于 完整 课题 
的 一 部 分 ,完整 课题 的 主要 研究 目的 是 探究 意大利 
蜜蜂 与 东方 塞 蜂 微 爷 子 虫 的 相互 作用 机 制 ,一 方面 
是 解析 宿主 的 胁迫 应 答 和 免疫 应 答 机 制 , 另 一 方面 
是 解析 病原 的 侵 染 机 制 。 鉴 于 东方 蜜蜂 微 雹 子 虫 的 
增殖 周期 尚 无 定论 ,以 及 处 于 增殖 过 程 的 病原 存在 
孢 原 质 、 裂 殖 体 和 成 熟 爷 子 等 各 种 状态 ,综合 考虑 后 
选取 东方 蜜蜂 微 抱 子 虫 感染 后 7 和 10 d 作为 本 研 
究 取 样 测序 和 数据 分 析 的 两 个 时 间 点 。 目 前 仍 缺 
乏 不 依赖 于 蜜蜂 宿主 的 东方 蜜蜂 微 孢 子 虫 体外 培养 
体系 。 本 研究 中 使 用 的 东方 蜜蜂 微 抱 子 虫 纯化 孢子 
是 从 患 微 孢 子 虫 病 的 意大利 蜜蜂 外 勤 蜂 中 肠 制备 而 
来 ,其 中 可 能 既 含 有 经 增殖 而 致死 宿主 的 孢子 ,也 仿 
有 未 侵入 宿主 细胞 的 孢子 。 鉴 于 目前 没有 区 分 此 二 
种 类 型 钨 子 的 技术 手段 ,我 们 参照 前 人 研究 中 采用 
的 接种 西方 蜜蜂 所 采用 的 孢子 浓度 和 接种 方法 
( Antúnez et al., 2009; Huang and Solter, 2013; 
Huang et al., 2016, 2018; Rodriguez-Garcia et al., 
2018) , fii E 3k LEE BRA ic HI Los MET Uf T 
虫 , 相 对 于 代谢 和 生命 活动 活跃 的 增殖 阶段 (如 裂 
殖 体 ) ,休眠 态 孢 子 中 仅 存在 较 低 水 平 的 代谢 和 生 
命 活动 (Williams et al., 2005; Corradi et al., 2008) 。 
本 研究 中 ,用 于 接种 工蜂 的 孢子 (其 中 已 完成 侵 染 
的 孢子 仍然 为 休眠 态 孢 子 ) 来 源 于 同一 批 制备 的 东 
方 蜜蜂 微 孢 子 虫 纯化 孢子 ,通过 对 照 组 和 处 理 组 数 
据 比 较 分 析 得 到 的 基因 差异 表达 信息 应 由 病原 侵 染 
所 致 。 本 研究 中 ,在 NeCK vs NeT1, NeCK vs NcT2 
FU NcT1 vs NcT2 中 分 别 鉴 定 出 1 397, 1 497 和 52 
个 显著 DEG ,表明 东方 蜜蜂 微 雹 子 虫 在 其 侵 当 意 大 
利 蜜蜂 工蜂 的 过 程 中 伴随 着 活路 的 基因 差异 表达 ， 
暗示 这 些 DEG 中 蕴含 着 病原 侵 染 相关 的 重要 信息 。 

Cornman 等 (2009 ) 组 法 及 注释 了 东方 蜜蜂 微 孢 
子 虫 的 基因 组 ,分 析 发 现 东 方 蜜 蜂 微 孢子 虫 的 供 能 
主要 依赖 碳水 化 合 物 代谢 、 糖 酵 解 、 氧 化 磷酸 化 和 磷 
酸 成 糖 途径 。 本 研究 发 现 ,有 10 个 基因 在 NeCK vs 
NcT1 , NeCK vs NcT2 £I NcTI vs NcT2 3 个 比较 组 均 
上 调 表达 ,其 中 4 个 上 调 基 因 涉及 碳水 化 合 物 代谢 ， 
包括 己 糖 激酶 基因 ( GenBank 登录 号 : 9423428) H 
TH -3 -BER Mi UE DS] ( GenBank 登录 号 : 9423035) 、 
磷酸 甘油 酸 激酶 基因 (CenBank 登录 号 : 9423279) 
和 丙酮 酸 脱氧 酶 el-B 亚 基 基因 (CenBank 登录 号 : 
9423243) (R 2) ,表明 东方 蜜蜂 微 人 钨 子 虫 在 侵 染 过 








































































































































































































程 中 通过 上 调 上 述 4 个 基因 的 表达 量 ,满足 自身 增 
殖 的 能 量 需 要 。 葛 麻 毒 素 是 一 种 可 抑制 蛋白 质 活 性 
和 改变 宿主 细胞 激素 编码 基因 的 表达 水 平 的 细胞 毒 
素 , 具 有 结合 并 粘 附 细胞 的 葛 麻 毒素 B 凝集 素 的 结 
构 (Spooner and Lord, 2015 )。 有 人 研究 表明 在 东方 蜜 
Me cf Y EER AE Bombyx mori Bl £& BmN 细胞 的 
WEEP, URE ZR B ERR BEES yn AR A pA fe. T E 
对 宿主 细胞 的 粘 附 能 力 (Liu et al., 2016b), Huang 
A (2016) 研究 发 现 东 方 蜜蜂 微 孢 子 虫 感染 西方 蜜 
蜂 的 次 日 就 能 检测 到 区 麻 毒素 编码 基因 的 表达 , 进 
一 步 分 析 发 现 葛 肪 毒素 含有 信和 号 肽 ,同时 宿主 的 核 
糖 体 蛋白 受到 抑制 。 本 研究 中 ,3 个 比较 组 共有 的 
上 调 基 因 包 含 3 PEREK B 凝集 素 编 码 基 因 
(GenBank 登录 号 分 别 为 9422649 9422645 和 
9422719)( 表 2) ,暗示 东方 蜜蜂 微 抱 子 虫 在 侵 染 意 
大 利 蜜蜂 工蜂 的 过 程 中 通过 增强 苑 床 毒 素 B 凝集 
素 的 分 泌 以 加 强 对 中 肠 上 皮 细 胞 的 粘 附 力 ,同时 掀 
制 宿主 细胞 蛋白 质 的 合成 。 细 胞 膜 蛋 白糖 基 化 是 细 
胞 识别 和 应 答 外 界 环 境 的 重要 方式 之 一 ( 刘 啸 尘 
等 , 2018) 。 本 研究 中 ,细胞 表面 糖 基 化 蛋白 编码 基 
(GenBank 登录 号 : 9423191) 在 NcCK vs NcTI , 
NcCK vs NcT2 £l NeT1 vs NeT2 3 个 比较 组 中 分 别 上 
381.74, 3.06 F 1.32 倍 ( 表 2) , HEM ART BEM 
子 虫 在 侵 染 过 程 中 通过 上 调 细胞 表面 糖 基 化 蛋白 基 
因 的 表达 对 宿主 细胞 内 环境 产生 应 答 。 此 外 ,共同 
上 调 基 因 中 包含 1 个 假定 蛋白 编码 基因 (GenBank 
登录 号 : 9422567) ,共同 下 调 的 1 个 基因 也 为 假定 
蛋白 编码 基因 (CenBank 登录 号 : 9422077) (R 2), 
说 明 东 方 塞 蜂 微 孢子 虫 的 基因 组 功能 注释 信息 尚 不 
完全 ,有 待 通过 分 子 克隆 、 基 因 和 和 蛋白 功能 研究 进 一 
步 完 善 。 到 目前 为 止 ,包括 东方 蜜蜂 微 孢 子 虫 在 内 
的 绝 大 多 数 微 抱 子 虫 的 基因 功能 未 知 , 最 大 的 制约 
因素 是 缺乏 转基因 操作 技术 平台 。Guo 等 (2016 ) 通 
过 碱 液 刺激 家 看 微 伯 子 虫 孢子 体外 发 芽 , 继 而 感染 
KÆ BmN 细胞 系 ,3 d 后 利用 脂 质 体 包 于 piggyBac 
转 座 子 载体 和 pIZT/ V5-His 非 转 座 子 载体 转 染 家 在 
微 孢 子 虫 感染 的 细胞 ,通过 连续 的 荧光 观察 和 分 子 
生物 学 检测 证 实 将 2 种 载体 携带 的 efp 基因 成 功 导 
人 家 在 微 抱 子 虫 基因 组 ;此 外 ,作者 还 通过 饲 喂 家 在 
微 抱 子 忠 孢子 感染 家 看 个 体 ,继而 注射 脂 质 体 包 吾 
的 上 述 两 种 载体 ,进而 也 通过 连续 的 区 光 观 察 和 分 
子 生物 学 检测 证 实 家 下 微 孢子 虫 基 因 组 插入 了 efp 
基因 。 目 前 ,笔者 团队 正在 参考 该 方法 尝试 建立 东 
方 蜜蜂 微 爷 子 虫 的 蔡 代 细胞 连续 感 当 系 和 转基因 操 
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作 技 术 体 系 。 

NcCK vs NcTI 和 NeCK vs NcT2 比较 组 包含 的 
DEG 富 集 在 代谢 进程 、 细 胞 进程 . 单 组 织 进 程 等 生 
物 学 进程 相关 功能 条 目 ,细胞 、 细 胞 组 件 和 细胞 器 等 
细胞 组 分 相关 功能 条 目 ,结合 .催化 活性 和 转运 活性 
等 分 子 功能 相关 功能 条 目 ,表明 东 方 蜜蜂 微 抱 子 虫 
通过 改变 新 陈 代谢 .细胞 生命 活动 相关 基因 的 表达 
水 平 促进 自身 增殖 。NeT1 vs NeT2 比较 组 包含 的 
DEG 主要 富 集 在 催化 活性 代谢 进程 和 细胞 进程 等 
功能 条 目 且 多 数 基因 上 调 表 达 , 说 明 东 方 蜜蜂 微 孢 
子 虫 在 侵 染 意大利 蜜蜂 工蜂 的 过 程 中 物质 和 能 量 代 
谢 活跃 。 进 一 步 对 DEG 进行 KEGG 通路 分 析 , 发 现 
NcCK vs NeT1 和 NeCK vs NcT2 比较 组 包含 的 DEG 
大 量 富 集 在 新 陈 代谢 相关 通路 ,包括 新 陈 代谢 途径 
(分 别 为 93 和 92 个 ) 次 级 代谢 产物 合成 (分 别 为 
37 和 39 个 ) .抗生素 合成 (分 别 为 29 和 31 A) IE 
叭 代谢 (分 别 为 21 和 21 个 )、 喀 啶 代谢 (分 别 为 20 
和 19 个 ) 碳 代谢 (分 别 为 1 和 18 个 ) 甘油 磷脂 
代谢 (分 别 为 11 和 13 A) 、 糖 醉 解 / 糖 异 生 (分 别 为 
11 和 12 个 ) 等 物质 代谢 通路 ,以 及 氧化 磷酸 化 (分 
别 为 10 和 10 个 ) 和 甲烷 代谢 (分 别 为 4 和 6 个 ) 等 
能 量 代谢 通路 ( 表 3)。 上 述 结果 再 次 表明 东方 蜜蜂 
微 孢 子 虫 侵 染 意大利 蜜蜂 工蜂 的 过 程 伴随 着 复杂 的 
物质 和 能 量 代谢 活动 ,病原 本 里 需要 提升 物质 和 能 
量 代 谢 满足 增殖 需要 ,但 同时 也 面临 宿主 的 抑制 。 
真菌 通过 MAPK cAMP 和 双 组 份 信号 通路 应 答 外 界 
环境 ,调控 其 生长 发 育 和 毒 力 REES, 2015), 
蜜蜂 球 圳 菌 是 男 一 种 常见 的 蜜蜂 真菌 病原 ,特异 性 

























































































































































































孢子 虫 的 侵 染 过 程 中 发 挥 重要 作用 ,下 一 步 拟 通过 
设计 相应 的 siRNA 对 GenBank 登录 号 分 别 为 
9424462 和 9422306 的 基因 进行 RNAi, 以 探究 其 
功能 。 

东方 蜜蜂 微 伯 子 虫 侵 染 蜜蜂 时 通过 中 空 的 极 管 
将 具有 侵 染 性 的 抱 原 质 注 入 宿主 的 中 肠 上 皮 细 胞 ， 
因此 极 管 蛋 白 在 东方 蜜蜂 微 抱 子 虫 的 侵 染 过 程 扮 演 
着 关键 角色 (Rodrfguez-Garcfa et al., 2018) 。 前 人 研 
究 报 道 干扰 极 管 蛋 白 3 编码 基因 的 表达 可 显著 降低 
东方 蜜蜂 微 孢 子 虫 孢子 的 数量 ( Rodriguez-Garcia et 
al., 2018) 。 本 研究 中 ,在 NecCK，NcT1 和 NcT2 中 
极 管 蛋白 1 编码 基因 (GenBank 登录 号 : 9423198 ) 
的 表达 量 分 别 为 171, 8 047 和 13 854, WERA 2 
编码 基因 (CenBank 登录 号 : 9423201) 的 表达 量 分 
别 为 333，17 919 和 27 591; 极 管 和 蛋白 编码 基因 
( GenBank 登录 号 : 9422340) 的 表达 量 分 别 为 190, 
2 923 和 3 994( 表 4) ,表明 东方 蜜蜂 微 雹 子 虫 在 侵 
染 增 殖 过 程 中 比 孢 子 状态 上 调 表 达 , 并 且 在 侵 染 过 
程 持 续 上 调 , 印 证 了 极 管 蛋白 对 东方 蜜蜂 微 和 抱 子 虫 
侵 染 的 重要 性 (Rodrfguez-Garcfa et al., 2018), $818 
BE APR AR Tr BE HE FH ELAR AS FR Hi SP ES 
良 环境 ,同时 在 应 答 外 界 环境 的 变化 传递 信号 和 附 
着 等 作用 (Southern et al., 2007), MWEE FA n] LJ 
CE EA HL FR TE HI , Bal Qe ec él fl) "b hS fen RE s 
AS 与 极 管 蛋白 2 和 极 管 蛋白 3 相互 作用 , 抱 壁 蛋 
白 5 对 极 管 蛋白 锚 定 在 孢 壁 上 起 到 非常 重要 的 作 
用 , 孢 壁 蛋白 9 定位 到 极 管 上 ,并 且 可 以 与 极 管 蛋 白 
1 和 极 管 蛋白 2 互相 作用 (Yang et al., 2018)。 本 研 

























































































侵 染 蜜蜂 幼虫 而 导致 白垩 病 ,笔者 团队 在 前 期 研究 
中 发 现 , 对 于 侵 染 意大利 蜜蜂 幼虫 的 蜜蜂 球 圳 菌 , 有 
48 个 DEG 富 集 在 MAPK 信号 通路 ,它们 的 表达 水 
平 随 蜜蜂 球 址 菌 胁 迫 时 间 的 延长 而 逐渐 增强 ( 陈 大 
福 等 ，2017 ) 。 本 研究 中 , 对 于 NeCK vs NeTI 和 
NcCK vs NcT2 比较 组 ,分别 有 7 和 7 个 DEG 富 集 在 
MAPK 信号 通路 ;笔者 前 期 研究 发 现 , 对 于 NcCK， 
NcT1 和 NeT2 ,分 别 有 5, 7 和 7 个 涉及 MAPK 信号 
通路 的 基因 高 量 表达 ( 熊 浴 玲 等 , 2019 ) ;深入 分 析 
发 现 本 研究 中 MAPK 信号 通路 上 有 2 个 DEG 与 前 
期 研究 的 高 表达 基因 重 码 ,其 中 1 个 DEG(GenBank 
登录 号 : 9424462) 在 NcCK, NeT1 和 NeT2 均 高 量 
表达 (FPKM fH 4) 3l] Jj 966, 829 Fil 22), 5 9h 1 S 
DEG( GenBank 登录 号 : 9422306) 在 NcT1 和 NcT2 
高 量 表达 (FPKM 值 分 别 为 209 和 205) ;上 述 结果 
表明 MAPK 信号 通路 及 其 富 集 DEG 在 东方 蜜蜂 微 




























































































究 中 ,一 共 检 测 到 5 个 孢 壁 蛋 白 编码 基因 (GenBank 
登录 号 分 别 为 9422782, 9422311, 9422437, 
9422653 和 9423533) .1 个 孢 壁 蛋白 前 体 编码 基 
(GenBank 登录 号 : 9422389 ) FI 1 AJE EE AN A AE E 
复合 蛋白 编码 基因 ( GenBank 登录 号 : 9424403) ,其 
HP JE BE TE A i 35 SE D] ( GenBank 登录 号 分 别 为 
9422653 和 9423533) 在 NeCK vs NcT1 和 NcCK vs 
NeT2 上 调 表达 ( 表 4) ,推测 它们 主要 在 孢子 侵 染 和 
增殖 过 程 中 起 到 重要 的 作用 ,而 孢 壁 和 锚 定 盘 复合 
蛋白 编码 基因 (GenBank 登录 号 : 9424403) 在 NcCK 
vs NeT1 和 NcCK vs NeT2 上 调 表达 ,推测 其 在 东方 
蜜蜂 微 孢 子 虫 侵 染 增殖 过 程 发 挥 重要 的 作用 。 此 
Db TEREA 9 编码 基因 在 NeCK 中 的 表达 量 很 高 
( FPKM 值 为 12 059) ,而 在 NcT1 和 NeT2 的 表达 量 
相对 较 低 (FPKM 值 分 别 为 2 807 和 2 885) ,暗示 其 
TEAR Ty E We a ed FH FSR R A pe 
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殊 作 用 。 本 研究 发 现 , 另 有 4 BREE ZR B 凝集 素 
编码 基因 (GenBank 登录 号 分 别 为 9422642, 
9423906, 9422644 和 9422647) 分 别 在 NeCK vs 
NeT1 ( 43923 4.77, 2.32, 2.72 $05.62 倍 ) , NcCK 
vs NcT2 (分 别 为 4.98, 3.23, 3.21 $403.38 倍 ) 中 上 
调 表达 ( 表 4) ,再 次 说 明 萝 床 毒 素 B 凝集 素 在 东方 
蜜蜂 微 孢 子 虫 侵 染 过 程 中 的 重要 性 。 

东方 塞 蜂 微 孢子 虫 寄 生 蜜 蜂 中 肠 上 皮 细 胞 ,由 
于 缺少 典型 的 线粒体 ,所 以 自身 只 能 通过 糖 酵 解 途 
径 产生 少量 ATP , 绝 大 多 数 的 ATP 必需 依赖 宿主 细 
胞 提供 (Cornman et al., 2009) 。 Huang 45 (2016) 通 
ib ZR; 38 epi fel-T- d I ES) Pg se T€ 1 -6 d 
H5 3B ET TUM FE , EMR 75 SEE PEL HE AE 
解 的 关键 酶 的 基因 在 感染 后 1 d 就 开始 表达 , 己 糖 
激酶 在 感染 后 2 d 被 检测 到 并 且 在 2 -6 d 的 表达 量 
上 调 , 糖 酵 解 途径 的 其 他 核心 酶 在 感染 后 3 d 也 被 
检测 出 来 ,同时 发 现 ATP/ADP 移 位 酶 基因 和 ABC 
转运 蛋白 基因 在 东方 蜜蜂 微 孢 子 虫 侵 染 过 程 中 也 维 
持 高 量 表达 。 本 研究 中 ,涉及 东方 蜜蜂 微 孢 子 虫 糖 
醇 解 的 几 种 关键 酶 基因 如 己 糖 激酶 基因 ( GenBank 
登录 号 : 9423428) \ 丙酮 酸 激酶 基因 (CenBank 登录 
号 : 9424331 ) 和 6- 磷 酸 果糖 激酶 基因 (CenBank 登 
RKF: 9424162) 在 侵 染 过 程 较 孢 子 状态 表达 量 上 调 
(#5) ,与 前 人 研究 结果 (Huang et al., 2016) 一 致 ， 
再 次 说 明 东 方 蜜蜂 微 移 子 虫 在 其 增殖 过 程 中 通过 糖 
酵 解 途径 提供 部 分 能 量 。 有 研究 表明 家 和 看 微 雹 子 虫 
Allie HH hf FHL Nosema locustae 的 已 糖 激酶 都 包含 
言 号 肽 , 均 为 分 泌 和 蛋白 (Timofeev et al., 2017; 陈 红 
丽 等 , 2018)。 本 研究 发 现 ,东方 蜜蜂 微 孢 子 虫 的 已 
糖 激 酶 同样 也 包含 信号 上 肽 ,推测 东方 蜜蜂 微 抱 子 忠 
在 侵 当 过程 中 将 己 糖 激酶 分 泌 到 宿主 细胞 质 以 调控 
宿主 的 糖 酵 解 。ATP/ADP fef Bi js i LT By HC 
宿主 能 量 的 关键 酶 之 一 ,Paldi 等 (2010 ) 通 过 RNAi 
对 东方 蜜蜂 微 爷 子 虫 的 ATP/ ADP 移 位 酶 基因 进行 
敲 减 ,发 现 东方 密 蜂 微 孢 子 虫 感染 的 西方 蜜蜂 体内 
的 病原 增殖 被 显著 抑制 。 本 研究 中 ,有 3 个 ATP/ 
ADP 移 位 酶 的 编码 基因 (CenBank 登录 号 分 别 为 
9424363, 9422880 和 9424586 ) 在 NcCK vs NcTI 
[ log, (Fold change) 值 分 别 为 0.72, 3. 60 $1 2. 31 | 
和 NcCK vs NcT2[ log, ( Fold change) 值 分 别 为 0.45 , 
3.44 和 2.34] 中 的 表达 水 平 上 调 ( 表 5) ,表明 ATP/ 
ADP 移 位 酶 对 于 东方 塞 蜂 微 雹 子 虫 客 取 和 宿主 细胞 
能 量 、 促 进 自身 增殖 的 重要 性 。ABC 转运 蛋白 的 主 
要 功能 是 结合 ATP, 利 用 ATP 水 解 产生 能 量 将 物质 










































































































































































道 浓度 梯度 运输 ( 冯 振 月 等 ,2018)。 本 研究 中 ， 
GenBank 登录 号 分 别 为 9422518 和 9423041 的 ABC 
转运 蛋白 编码 基因 在 NeCK vs NeT1 [ log, ( Fold 
change) 值 分 别 为 0.93 和 1. 83 ] 和 NcCK vs NcT2 
[ log, (Fold change) 值 分 别 为 0.64 和 1. 85] 中 上 调 
表达 ,体现 出 对 东方 蜜蜂 微 抱 子 虫 侵 染 过 程 中 物质 
运输 的 重要 性 ;但 另外 4 个 GenBank 登录 号 分 别 为 
9423050, 9423588 , 9422298 和 9422297 的 ABC 转运 
蛋白 编码 基因 在 NcCK vs NeT1 [ log, (Fold change) 
(3g -3.65, -1.72, - 1.91 和 -0.19] 和 
NeCK vs NcT2| log, (Fold change) fH 13] Jy — 3.34, 
-0.13, -1.56 和 -0.35] 中 的 表达 水 平 也 为 下 调 
( 表 5) ,鉴于 东方 蜜蜂 微 雹 子 虫 与 蜜蜂 之 间 存 在 密 
切 的 相互 作用 ,上 述 3 个 基因 的 下 调 表达 可 能 是 由 
于 东方 蜜蜂 微 和 抱 子 虫 受到 宿主 的 抑制 ,具体 的 机 制 
有 待 于 进一步 深入 研究 。 

AS RIED SAR Tr BEE HLF HEALY EKA 
蜜蜂 工蜂 过 程 的 差异 基因 表达 谱 , 在 mRNA 组 学 层 
面 揭示 了 东方 密 蜂 微 孢 子 虫 的 侵 染 机 制 ,筛选 出 病 
原 的 毒 力 因 子 和 侵 染 相关 因子 基因 可 为 后 续 功 能 研 
究 提供 候选 靶 标 ,并 为 治疗 蜜蜂 微 雹 子 虫 病 提 供 湾 
在 的 分 子 靶 点 。 
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